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Abstract 

By suitably re-scaling the conformal Einstein's equations we are able to apply recent 
results in the theory of PDE, and prove that they possess slow solutions in a future 
neighborhood of an initial surface reaching X + . The structure of the equations obtained 
allows to split (up to any given order) the initial data into those generating slow 
solutions, i.e., those driven by the sources, and those generating fast solutions, i.e., 
those which represent gravitational radiation with no relation to the sources. Thus 
effectively resulting in a proposal to prescribe initial data for solutions with no extra 
radiation up to the order needed for each given application. 



1 Introduction 

A very important step in the study of generation of gravitational waves by the self gravitating 
motion of celestial bodies or black holes is to determine which initial data must be given in 
order to start the computation of the gravitational dynamics. For one would like to specify 
data with the minimal amount of gravitational radiation not related to the sources, that is to 
tailor the initial data for gravity to the data for the sources. This task is particularly difficult 
in general relativity for the nonlinearities of the theory do not allow for a clear cut resolution 
into source driven and source independent solutions. Fortunately a theory describing in 
detail the behaviour of solutions to equations with different time scales has been developed 
in the last two decades which allows to approach the problem in a more systematic way. For 
the case of the problem at hand, which we have taken to be an asymptotically flat space-time 
with matter sources on it, there are two clearly distinct time scales, a slow one which can be 
associated to the characteristic Newtonian time scales, for instance the orbital period of a 

* Supported by STINT, The Swedish Foundation for International Cooperation in Research and Higher 
Education. 

tFellow of Se.CyT-UNC. 
tMember of CONICET. 



1 



test particle circling around the sources near their boundaries, that is a distance associated 
with its mass and size. The most important characteristic of this time scale is that does 
not change if one increases the value of the speed of light. The other time scale, faster than 
the previous, is the time light takes to travel across the source, this clearly changes as the 
speed of light is increased. Thus, by changing the value of the speed of light in the equations 
one is able to distinguish between these two time scales. In turn these two time scales can 
be used to distinguish between solutions carrying gravitational radiation generated by the 
sources from those which represent pure radiation, not coming from any source. The former 
would move without appreciable change as we take the speed of light to infinity, while the 
latter would oscillate faster and faster, thus in the limit becoming singular. But does this 
difference between solutions exists? and how can we characterize it? Physically we know 
that something like this must happens for this is what we impose when we speak of the 
Newtonian limit, and subsequent postnewtonian approximations, and we know they are a 
good approximation to describe some celestial systems. Mathematically the existence and 
characteristics of this splitting of solutions, which is only asserted up to any finite given 
order, is provided by the theory of different time scales mention above, see for example 
i and §. 

For these mathematical theory to apply it it must be assumed that the system of partial 
differential equations ruling the dynamics is symmetric hyperbolic and become singular in 
a very specific way. If we define the parameter controlling the ratio between the two time 
scales as e, in our case the inverse of the speed of light, the systems must be of the following 
form 

(1.1) A ° i3 (eu k )^ = + A a li3 {u\e))d a u> + B^u^e), 

with the matrices K a A± and A symmetric, being K a constant, A positive definite and the 
matrices A\ and the vector B continuous in e and uniformly continous in u for bounded u 
(see 1, [!,!). 

If the system is as above then one can assert that there are one parameter families (in e) 
of initial data sets which, for a finite period of time, give raise to solutions whose derivatives 
up to some given previously specified order are bounded in the limit e — > 0. Those are 
the slow solutions one is seeking in order to study the gravitational radiation produced by 
sources. Even more, the theory also tells how the rest of the solutions, called fast solutions 
would behave when they are small (order e) compared with slow solutions. The choice of 
initial data giving slow solutions is called initialization and, as the theory shows, simply 
consists in making sure that the time derivatives of the solution up to the desired order are 
bounded at the initial surface. Using the evolution equations these requirements become 
conditions on the initial data and its spatial derivatives. 

In [|7| and || the task of writing Einstein's equations in a form like |(1.1)| was carried out 
in the context of a time function associated with an asymptotically flat foliation of space 
time, thus the assertion that the solutions arising from those special initial data sets had 
regular limit did not hold all the way up to null infinity, and so the results could not be 
applied to obtain conclusive estimates about the gravitational radiation produced. 

In order to truly study the radiation phenomena, in this paper we rewrite Einstein's 
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equations in an asymptotically null foliation of space-time, and accomplish the task of casting 
them in the form |( 1 . 1 )| . This is done in the next section, where we heavily use Friedrich's 
schema of conformally regularizing Einstein's equations (c.f 0). Once this is done we use, 
in the third section, the resulting system of equations to obtain the Newtonian limit, which 
in this context tells how to select initial data sets whose solutions are bounded for a finite 
time intervall, that is we perform a first order initialization of the data for the system. It is 
worth noticing at this point that besides evolution equations, the Einstein system has also 
constraint equations between the initial data. They complement and are consistent with the 
ones the initialization scheme produces, leading together into the Newtonian limit equations. 

2 The conformal Einstein field equations as a symmet- 
ric hyperbolic system 

In order to dimensionalize the Einstein equations we introduce, as in |7|], the concept of 
dimension for geometrical objects. We impose that the coordinate functions have dimension 
L and define the dimension of the components of a vector n a , [n a ], such that when the vector 
acts on any function /, it gives a function of dimension [n(f)) = [n a ]j-[f]. We define the 
dimension of a covector m a , [m a ], such that when the covector acts on any vector n a it gives 
a function of dimension [m a n a ]. We extend these definitions to tensor fields in the obvious 
way. This implies that the metric tensor is dimensionless, the connection has dimension of 
L^ 1 and the Riemann tensor has dimension of L~ 2 . The parameter which shall distinguish 
between slow and fast solutions is, e = c _1 , where c is the velocity of light. 

We shall deal with a four-dimensional Lorentzian manifold (M, g a b)f\, where g a b has sig- 
nature (— , +, +, +) and satisfies the Einstein's field equations 

(2.1) R a b — -Rfjab = KE 4 T ab 

with k = 8ttG, G been the gravitational constant. Since the dimension of the Ricci tensor 
is L~ 2 , the energy-momentum tensor T ab has dimensions of energy density. 

Furthermore the source is thought to have compact support, physically we are thinking 
in a isolated system, for instance a binary system. Thus we shall deal with asymptotically 
flat space-times at future null infinity fllT| |, i- e - a triple (M,g,Q) that satisfies: 



1. M is a four- dimensional manifold with boundary X + , 1 + is diffeomorphic to R x S 2 . 

2. g a t is a Lorentz metric on M, (M, g a b) is time and space oriented and strongly causal. 

3. Q is a function ('the conformal factor ') on M with 

Q > on M = M\J+; Q = 0, dQ ^ on J + . 

4. T + is a null hypersurface with respect to g a b in the past of M. 



1 Here, we use the abstract index notation, see Col 
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5. The metric g ab = 2 gab satisfies the Einstein's field equations 



6. There is a spacelike hypersurface 5 of (M,g) with two-dimensional boundary Z which 
also belongs to X + . The hypersurface S is diffeomorphic to the closed unit ball in R 3 
whence Z is diffeomorphic to the sphere [] S 2 . 

Expressing the field equation |(2.1)| for g ab in terms of g a b and Q, and splitting them into 
the trace-free part and the trace yield 

1 K 4 o 

(2.2) Qa ab = -v a v b n + -g ab v c v c n+'^-nTab 

(2.3) n 2 R = -K£ 4 T-6(fiV c V c fi-2V c fiV c fi), 

where R is the Ricci scalar, la ab the trace-free part of the Ricci tensor and V denotes the 

o 

Levi-Civita covariant derivative with respect to g ab , Tab is the trace free part of the energy 
momentum tensor field and T = g ab T ab . Remark that since the sources are assumed compact 
support, T ab has a natural extension to M. 

Under the rescaling (g,Q) — > (G 2 g, Gf2), where the function Q is positive everywhere, 
the Ricci scalars of g^ v and of Q 2 g^ u are related by 

e R[g^} - G 3 R[Q 2 g^) = 6 v A v A e. 

Thus given any Cauchy data for G on an initial surface S with G positive on S, one may 
determine G in a neighbourhood of S such that the Ricci scalar of the rescaled metric 
becomes constant. We shall use this conformal freedom to set the value of R to be 6s 2 . 

As a first step to treat the conformal Einstein field equations as a system with two different 
time scales, we cast them as a first order symmetric hyperbolic system with constraints, using 
the formalism developed by Helmut Friedrich in |J. 



2.1 The Choice of Frame 

Following Friedrich, we choose an orthonormal frame, this choice yields an hyperbolic system 



of the form (1.1) , where the matrices K a depend on the frame. Since we want to describe 
the Newtonian limit, i.e. the limit when e — > 0, we need that the matrices K a be constant 
in that limit on the hypersurface S on which we give initial data . We achieve this result 
demanding that the frame become constant in this limit. 

To give a better idea of this choice of frame we exemplify it in Minkowski space-time. 
The line element takes the form 

dS 2 = -^! + dr' 2 + r'W, 
e 2 



2 Most of our resuls are local, so some topological conditions are superfluous. 
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where du 2 denotes the standard line element on the two-dimensional unit sphere. A possible 
hypersurface on which we could give initial data is the spacelike hyperboloid given by 

S = {t' 2 - e 2 r' 2 = 1, t'>0}. 

By performing the coordinate transformation 

't + er\ , , ft — er 



t' + er' = tan (^-^-) , t' ~ er ' = tan (- 



and rescaling the line element with the conformal factor 

r^, \ c ft + er\ ft — er 
\l{t,r) = 2 cos I — - — I cos 

one obtains the unphysical line element 



(2.4) ds 2 = -^ + dr 2 + S ^f±du 2 

e l e 2 



Minkowski space-time is only the region 

< t < 7r, < r, \t + er\ < 7r, \t — er\ < it. 

In this new coordinate system the closure of the initial surface S is given by S = {t = 
7r/2, < r < 7r/(2e)}, this surface represents a ball of radius 7r/(2e) and when £ — > it 
becomes the three-dimensional Euclidean space. Note that the surfaces t' ± er' = const. 
represent the advanced and retarded null cones and in the limit (when e — > 0) they become 
the hypersurface t = const., where t represents now the absolute time. 

We claim that we can choose on S coordinates such that the frame becomes constant 
when e — > 0. To see this we note that in Cartesian coordinates the induced metric on S is 
of the form 

ds 2 = dx 2 + dy 2 + dz 2 + e 2 h ij (e, x i )dx i dx j , 
for a known smooth in (e,x l ) matrix hij. Thus the frame to be choose is of the form 

d , d 



eo = £tt7> ei = e 3 - 



df 1 l 8xi' 

where c l j = 5 l j + ee l j, for some smooth function e % j. 

Notice that the second fundamental form associated with this surface, Xiji vanishes, the 
conformal factor is Vt(n/2,r) = 2 cos ^ 7r / 2 + £r ^ cos ( j^ 2 ~ er ^j and defining S = e (Q) on S we 
have S = —e. This complete the example. 

We now turn into the general case, it is clear that for metrics close to Minkowski we can 
choose an orthogonal frame e a ^, with /i — 0, 1,2,3, in such a way that the components of 
the metric in this frame become 
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Given a spacelike hypersurface S, we choose eo as the unit normal to the hypersurface and 
three orthonormal vector fields (i = 1, 2, 3) on S, which will be propagated parallel in the 
direction of eo- Furthermore, we fix Gaussian normal coordinates x a on a neighborhood of 
S in D(S), hence 



where are tangent to the hypersurfaces S t — {t — const.} and e J i| £= o = 5^ 

We denote by 7^ the connection coefficients of the canonical covariant derivative V 
with respect to e a V) i.e. 

V e a = ^ e a 

and with the choice of the frame given above, the connection coefficients have the following 
properties 

7o/ = 0, 1^ V 9\) V = 0. 
Finally we choose the gauge source as in the Minkowski case, i.e. R = Qe 2 . 



2.2 The Variables 

Following Friedrich's work and after several steps which are detailed in the Appendix, the 
conformal Einstein equations can be split into constraint and evolution equations for the 
variables 

u (e Q,e j, 7 iji Xiji ^> ^0? @oii Eiji -Bjj), 

where QEij and QBij represent the electric and the magnetic parts of Weyl tensor, Xij — 
represents the second fundamental form on S. The variables So,Sj and s, when the conformal 
Einstein's equations are satisfied, become T, v = V„f2 and s = ^V^V^fi respectively. Note 
that the variables have been already rescaled in such a way that the source of the system 
becomes regular in e and the Newtonian potential appears in the first initialization. Thus 
these variables correspond to the "tilde" variables in the Appendix. 



2.3 The Conformal Evolution Equations: 



n t - s = 

+ ^ ( i SEo + J% + J_s ( ) _ i Eo = o 

3 3k£ 2 ke a 

Eo )t + e 2 fia 00 + + s 5— - "^^^ = 

4 8 8 

3 

£ M + e 2 VLa 0i + ^-Q J, = 
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e\,t + X l i + ex 3 1 e l 3 = 

(2.5) e\ t = 

Xij,t + £XijX l i ~ eQEij - e (<Jij - gijcroo) = 

fjk,t + £ 2 YikX l j - eVLB jT e Ti k + e (#Vo* - 9kj^o) = 

&Oi,t - ~V 3 <Tij = 

e 

&ij,t — -V(i<7j) — Z Eij — eTii B T (i e Tl j) — -VLKt (ij) = 

Eij,t + -DkBi(i€j) kl — 3eE k (iXj)k + ^Xk k Eij + egijE km x km + ~«*o(»j) = 

1 k 

Bij,t D k Ei(i6j) — 3eB (iXj)k + %£Xk B {j + eg^B^x m + —e n (jtj) n i = 

where D k is the covariant derivative on the hypersurfaces S and p, J 1 , E 2 S l i and t^x, 
represent, respectively, the energy density, the momentum density, the pressure tensor and 
derivatives of the energy momentum tensor, as seen by an observer at rest with respect to 
the foliation chosen. They are defined by 

(2.6) T ab = p e\e\ + 2e J {a e b) + e 2 S ab , 
and 

(2.7) txrjv := ftT 1 ^ [\F V ] V = f2 _1 (V[\F V ] U + -1 E[ A i^] w — g^F^pE^) , 

being F vu = | — ^g vu T^j. In spite of the presence of | in the source terms of the three 
last equations of |(2.5)| it can be seen that they are smooth in e (see Appendix). 

This system is of the form as can be seen after a lengthy and mostly uninteresting 

calculation. Thus if appropriately coupled with symmetric hyperbolic equations describing 
the matter fields 3 it results in a bigger system with the same properties, see for example |]12" 



Thus the general theory of systems with different time scales can be applied in a straight 
forward way. 



2.4 The Conformal Constraint Equations: 

The projection of the conformal Einstein equations on S yields the following constraint 
equations: 

DiQ - e 2 ^i = 

hhS^ K/£^ o KiS^ t\j&^ 3<E^ 

DiS - £ 3 o- i0 £ + eVyE 3 Y JiT>0 2~ Sij ^ ~~ 24tt DiP + 24tt DiS + = ° 

3 Here we are assuming two things, first that the matter equations are not singular in e and second 
that boundary conditions for the matter fields can be handled appropriately in order for the usual energy 
estimates to hold, (the case of electromagnetism can be treated without further problems by adding it to 
the singular sector). 
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A.So - e 3 Xij ^ + e 3 n(T i0 + —^IJi = 

KE 4 o 

D.-X'j ~ £Xij^o + e^a,^ + —gij - sg i:j — Q Sij= 

(2.8) 2e\ kJ] + e^/j - e)/ k ) - e(f jk - 7 \y) = 

3 r l jM + 2e 2 x\kXi]j - 2e 2 Q{8\ k Ei]j - lj\ k Ei\] - 2e 2 5\ k ai] j - 2e 2 gj[i<j l k ] - 5\ k gi}j = 

2-D[jXfc]i — s 2 VL B iT e T jk — 2e 2 gi^a k ^ = 

2D[i(Tj]o — 2E 2 x k [i(Jj]k — £ 2 £fc B k T e T ij — Vlnt^ = 

2D[i(Tj] k — 2e 2 Xk[i^j]o — 2£ 2 ^i{5 l [iEj} k — l k ^Ejj 1 } — £:S B kr e T \j — VtKt ijk = 
D l Eij — e 2 B ik x % i £ hl j + k hjo = 
D l Bij + e E ik x\ e j + ~z t mn0 e mn j = 0, 

where 3 r* j k \ is the Riemann tensor of the hypersurface t = const, given by 

Vj-fcj = e k {fij) - ei(f kj ) + e 2 (^ km l m ij - l\ml m kj - Y mj (l m ki - l m ik)) 

o 

and Sij is the traceless part of SV,-. We shall assume for the rest of the paper that there exist 
smooth in e one parameter families of solutions to the constraint equation that are C p , for 
some p > 1. The existence and regularity of the solutions to the constraint equations have 
been studied in [0, and the smoothness of the solutions in the parameter e should follow 
from the implicit function theorem. 



3 The Newtonian limit 

In the present context the Newtonian limit is defined as the subset of the one parameter 
family (in e) of solutions to the equation system |(2.5)| - K2T8)| which are bounded as e goes to 
zero. 

According to the theory of different time scales, to achieve this result it is enough to 
initialize the equation system, i.e. to require that the time derivative of the dynamical 
variables on S remains finite when e goes to zero, being this requirement consistent with 
the constraint equations. Since the solutions will remain bounded in time, as the theory of 
different time scales asserts, it implies that these solutions will satisfy similar initialization 
and constraint equations on each future hypersurface of the foliation. 

Assuming the dynamical variables are C 1 in e, and that the initial data can be set as 

u(x,0) = u {x) + 0(e), 

the singular part of the evolution equations |(2.5)| with the above requirement yields at the 
initial surface 4 S 

D j aij = 



4 Since these equations will be valid for each slice S t , t € [to, T]. 
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(3.i; 



D l cr j = 
DkEi(i€ kl j) = 







and to zeroth order the constraint equations (2.8) become 



DiQ = 
D iS = 

DiEj = (s - -)g iS 
AS = 

(3.2) e [hJ] = 

3rl jki = ^ l [k9i]j (Gauss' equation) 
DyXk]i — (Codacci's equation) 
D[i<Tj]o = 

D[i(cr j]k - 9j]k^p) = 
D\E, l3 - ^g ijP ) = 
D l B tj = 0. 

Notice that the e = limit of the geometry at which we solve the initialization equations 
is rather singular. If we picture the compactification as immersed in the Einstein Universe^], 
then the part of the Cauchy surface of the Einstein Universe corresponding to our initial 
surface, becomes bigger and bigger, and in fact each Cauchy surface of Einstein Universe 
with topology 5* 3 in the limit blows up to become R 3 . 

We stress that the general relativistic dynamics, which is what we are interested in, after 
all, has a non-vanishing e. Thus, we shall be interested in the solutions to the initialization 
equations in just a portion of R 3 , that part that corresponds with our initial surface. 

If no boundary conditions are imposed on the solutions to the initialization equations an 
infinite set of solutions to them would be possible and so it would result in a non unique 
Newtonian limit. We believe this undeterminancy is due to the fact that the initial sur- 
face is not a Cauchy surface and that most of these solutions would in fact not arise from 
solutions which are asymptotically flat at space-like infinity, for they would not have the cor- 
rect (mild) singular structure there. We conjucture that the only solutions arising from an 
asymptotically flat space at space-like infinity are those which when considered as solutions 
in the whole limiting initial slice (R 3 with a flat metric) have asymptotically flat boundary 
conditions and, in what follows we shall treat just that case. 

5 This is just a convenient picture and is not totally correct, for in the presence of matter we know that 
at least the point corresponding to i° must be singular 
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We solve now the system |(3.1)| - K3~2~)| . The conformal freedom in the initial data is fixed 
as in Minkowski space-time, namely 

3 

(3.3) Q = 1, s = —-, and £j = — Xj. 

When we fixed the value of the Ricci scalar we still had a freedom left in the choice of the 
time derivative of at the initial surface S and since the mean curvature of S changes as 



dx 



o 



under a conformal transformation, we choose 7: — , and in turn E , in such a way that the 

rescaled mean curvature becomes zero. 

Consider now the equations for Xij, fixing j, these equations say that Xij = DjYj. Because 
of the symmetry of Xij we have DuVj] = which in turns implies that Vi = Di ip, for some 
arbitrary function ip. Thus if x — 0, then 

Aip = 0, 

imposing the asymptotic condition that the field vanishes at infinity, ip must be zero. 

Since Di<j j = 0, imposing the same asymptotic condition as above, we must have aoj = 0. 

Finally, under appropiate boundary conditions we claim that the unique solutions of the 
equations for B^, aij and Eij in |(3.1)| and |(3.2)| are given by 

Ik 1 

B^ = 0, = --DiDj(f) a + Qij-^P and = DiDj(p E - -g {j A(f) E 

where the functions (fi a and (pE satisfy 

A0 CT = 2nGp. 

and 

A(f) E = 2nGp, 

respectively. Uniqueness of the solution in the Laplacian equation yields <fi = <pE — 'Pa- 
We call this function the Newtonian Potential, and the initial data we have determined the 
Newtonian data. 

It remains to prove the uniqueness claimed above for the initialization, hence of the 
Newtonian limit. We do that by showing that the solution for a^, and B^ are unique. 

It is easily seen that the difference between two different solutions for the equation for 
is a symmetric tensor that satisfies D % 5aij = and D^Sa^k = 0. As for Xij, the symmetry 
of and the last equation imply that daij = DiDjip. The first equation implies that 

ADiip = 0, 

imposing the asymptotic conditions we have that if) = const, and so that 5aij = 0. The 
uniqueness for E^ and B^ follows from the next lemma as applied to the difference between 
two solutions 
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Lemma 3.1 Let Vij be a smooth tensor field in R 3 such that Vu — > quickly enough when 



oo. Moreover let Vy be a symmetric, transverse and traceless tensor satisfying 



(3.4) 



D k V Ki e kl j} = 0. 



Then Vy = 0. 



Proof: Using de Rham's Theorem we note that if is a transverse tensor (i.e its divergence 
is zero) then there exists a tensor field W nm such that 



(3.5) 



D k W H e kl r 



Since is symmetric, |(3.5)| implies 

(3.6) Di W nl -D n W = 0, 
where W = g^ Wij. Inserting |(3.5)| in |(3.4)| we obtain 

(3.7) A W(ij) - D m D (t W Hj) = 0. 

Since and Wij satisfy linear equations it is sufficient to prove that Vij(W s ) = V{j(W ) = 0, 
where W s (W A ) is the symmetric (antisymmetric) part of W. 
We consider first Vij(W ), equation |(3.5)| implies that 







where k 



Imn pr/ 



'Di kj. 



e— W A mn . Then 

V^(W A ) 

The symmetry of V^j implies that kj = Djip for some scalar function ip, and since Di k l = 0, 
ijj satisfies A^ = 0. Thus ip = const. , which in turns implies Vij(W A ) = 0. 

Consider now the symmetric part W s . Using equation |(3.6)| , equation |(3.7)| becomes 



(3.8) 



A Wfj - Di Dj 



0. 



Let 4> such that Acf> = W s , where W is the trace of W$, then |(3.8)| becomes 



A(Wg- DJ) 



using the asymptotic conditions, uniqueness of solutions to Laplace equation implies, 



DiDj( 



Inserting this expression for W s into [[33], we obtain Vij(W s ) = 0. □ 
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4 Conclusions 

By suitably rescaling Friedrich version of Einstein equations we have cast them in a form 
where standard results of the theory of different time scales apply. In particular this theory 
asserts that there exist slow solutions, that is solutions which move at the pace of the matter 
motion and so describe the gravitational radiation they produce. Even more the theory tells 
how to single out these solutions by imposing conditions the initial data sets must satisfy. 
These conditions are usually called initialization conditions and they are naturally consistent 
with the constraint equations of General Relativity. 

Imposing these conditions to first order we find what in the present setting should be 
called the Newtonian limit: all degrees of freedom are frozen, only remains a scalar function, 
the Newtonian Potential, which is completely determined by the matter content of the space- 
time. 

These initialized solutions remain slow all the way into a finite neighborhood of T + , thus 
the radiation they imprint at X + can be truly assigned to the slow motion of the matter 
sources. 

The structure of the equations used allows to control the slowness of the solution to any 
order on e, by continuing the initialization procedure i.e. making sure that further time 
derivatives of the solution at the initial hypersurface remain bounded. These postnewtonian 
orders should give enough information as to make a rigorous justification of the Quadrupolar 
formula, in a similar scheme to the one used by Winicour on null hypersurfaces ||14|| . But 
we believe that the most important application to the above system of equations should be 
in numerical calculations for there, besides giving a symmetric hyperbolic set of equations, 
a result already know from Friedrich's work, one would have a tight control on the initial 
data to be prescribed. 



5 Appendix: The dimensionalized field equations 

In order to have the Conformal Einstein field equations as a symmetric hyperbolic system 
with smooth sources in e we shall rescale the variables used in [Pi. These variables are 



u 



-j a 



[11 T' vr\ i ^Jj/, S, (J^v , d^urjX ) j 



where Vtd^ ur) \ is the Weyl tensor. In these variables the equations |(2.2)| and |(2.3)| becomes 



(5.1) 
where 



v\r\) 



M = V M fi-S M 



KE A 



T a vr, Y 



KE 
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(5.2) Qfxv = V^S;, + {~lo~^ w — sg^y — Q Tpv 

T" vX <* = (7^-7Me;-(el/-e:/) 

E^ v Xr\ ^ v\r\ R^uXtj 

L„\ v = 2V[\cr ri ] v — Y>^ v \ v — VtKeH\ w + — V[xRg v ]v 
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where R^ v xr) denotes the Riemann tensor defined by the metric, and r fJ- u x ri the Riemann 
tensor defined by the connection and the torsion T^ uX , i.e 



Besides the unknown u, the function R appears in the equation system z = 0. This is a 
gauge source function for the conformal factor which shall be chosen in what follows to be 
constant as in the Minkowski case, i.e. R = 6e 2 . 

Qui* = entails that the trace free part of the conformal Einstein's equation (2.2) is 
satisfied, T^ v \ = and K^ v xr, = guarantee that the metric reconstructed from the frame 
is the canonical one. L u x v = and H v \ v = are the Bianchi identities with matter and 
Pp = is an integrability condition for Q^ v = 0, i.e. it is obtained by taking the covariant 
derivative of Qp„, contracting and using L\ w = 0(see fL2| ). 

Instead of the variables we shall use the variables E ra and B ra . These tensor fields 
can be thought as the electric and magnetic parts of the Weyl tensor, i.e. 

(5.3) Vt 1 C flv x p = 2{l fl [xE p ] u — IvixEp]^ — n[xB p ] T e T p U — n\pB v ] T £ T \ p } 

where 



E T<T 


= Sl- l C^x P h\n v h\n p 




= n- l C; vXp h\n»h\n? 


I /IV 




hp,u 


= 9pu + n^n u 






pv\p 


— 2 ^pvafi e \p 



here we shall take n p to be the unit one form normal to the hypersurfaces t = const., i.e. 



n a = e a . 



Equation |(5.1)| can easily be split into the propagation equations along e and into the 



constraint equations projected on S. 

The Evolution Conformal Equations: By the gauge assumptions, Vo = £dt, so 
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= O = efii-S 

4 5 2 

= P = es, t + a^gT ~ ^ T o„ + + y^o 

= = £^i/,t + fio"o^ — ^fiW ^~ ^ To* 

= T% x e; = eel t + Yxoe; 

(5.4) = X^o/3 = £7V,t + 7V7%-^(-eo^ + eo^/ + J B / 3re r ^)- 

-25 M [0 o-/3]i, - 2^[^o- M ] - £ 2 ^[o 
= -g vX L vXi ^e(T^ t -V j a ij + n^g v h Xiv 
= i^joj = £Vij,t — VjfToj — TjqEij — YnBi T e Tl j — QneHoij 
= #(i|o|j) = + D k Bi(iej) kl - 3E k (iXj)k + ^Xk^ij + 

+gijEkmX hm + KeH (ij) 

= —-H n i^e nl j) = eB ijtt — D k Ei^e,j^ 1 — 3B h ^Xj) k + 2xk k Bij + 

+gijBkmX m + —£ £ n (itj)nl- 

The Conformal Constraint Equations: 

= Oi = DiQ - Si 

= Pi = DiS + a iu ZP - — Tiu + ^—DiT + — E; 

2 24fi 2 

= Qio = A S - + Qa i0 — Tj 

= Q^- = DiHj — Xij^o + Q&ij — sgij — £ 2 T^- 

= T* ifc = ( 7 V - 7%) - (el/j ~ e)/ k ) 

(5.5) = K i jkl = 3 r i jkl + 2 X \ kXl]j -2n{5 i [k E l]j -l j[k E l] i } 

-2S\ k anj - 2g jV (j l k ] - e 2 5\ k g t]j (Gauss' equation) 
= K° ijk = 2D{jXk}i - QB iT e T jk - 2g i { j a k } 

(Codacci's equation) 
= L 0i j = 2D[ i aj] — 2x k [iOrj] k — Y lk B k T e T i j — {lKeH i j 
= L kij = 2D[i<jj] k - 2xk{i&j]o - 2T, lx {h' J '[iE j ] k - l^Ejf} + 

— S B kT e 7 ij — VLke tij k 
= H 0mn e mn i = D l Eij — B ik x\e kl j + K,e A t j 

4 

= H j = D l Bij + E ik x' l ie j + ——e mn jt mn0 , 
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Since the variables Bij and Eij in this gauge (c.f. ]T5|) satisfy a symmetric hyperbolic 
system of differential equations, the system above is, as in the case with the variables d l jki, 
a symmetric hyperbolic system with constraints. 

The sources of this system are not smooth in e, therefore in order to apply the existence 
theorems, we rescale the variables as follows. 

~ _ S o_o~- v ^° ~i 

s — ~i " — 9fiu — 9[Wi ^0 — i — ~i 7 jk — 



°V ~ 4 + £ 2 ' 



Part of this rescaling is based on the behaviour of some of these variables in Minkowski 
space-time with metric given by (2.4) . With it the Conformal Evolution equations become 



= fit -S 

KS 2 ( ~ 1 

= s it - £ 2 a 00 S + e 3 a 0i Ejg v + — (^3pS - ^p >t ) H 

" 1 n -st + J% + -Ls, t ) - % 



2 V4 12Q '7 4 

= S , + e 2 Qa 00 + ^ + s - ^V x p - ^V 1 ^ 

4 o o 

= S M + e 2 VLa m + ^-ttJi 

= e*i, t + xS + ej^ e*i 

(5.7) = eV 

= Xy,< + - efi-Ey - £ (ofy - ftj5- 00 ) 

= 7V, t + £ 2 fifcx'i - eVLB jT t Ti k + £ (^orofc - g k j^o 

e e 

1 _ _, _, _, { i _, 

= Oij.t — ~V(i(Tj)o — So-Ejj — eYji B T (i e T j) — -fi/ttofe) 

= £t i>t + -DkBuiEjf 1 - 3eE k {i Xj) k + 2ex fc *£y + E~gijE km y* m + |/cto«i) 

1 _ _, K 

= Si 3 - |t — -D k Ei^ej) — 3eB (iXj)k + 2£Xfc -Bjj + egijB km x m + — e™ (ji,>« 

and the Conformal Constraint equations, 
= DiQ-e 2 ti 

4 6 2 4 q 2 

= DiS- e 3 a i0 t + e%jt? - —Jit - — <L S J ' - -^Ap + ^DiS + — S; 
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(5.8) 









4 

KB ~ 



2 ~ _ VI _ KE 4 ~ o 

DiHj - exij^o + e \la i:j + —gij - —il S 



'j 



2e 



-2~i ~ 



si d 



2D[jXk]i — Ay £ T jk — 2£ 2 gi\j&k]o 

2D[i(Tj]o — 2s 2 x k [id'j}k — £ 2 E k B k T e T ij — f2«t^ 

2D[iaj] k — 2e 2 Xk[iO'j]o — 2e 2 Tii{8\iE^ k — l k ^E^} — eS B kr e T ij — VLn,tij k 



S\k9i]j 



D'Eij — e 2 B ik x\ c l j + k tojo 



D l Bij + e 2 E ik X\ C kl j + — tmnQt 



K 



J- 



We claim that in the system |(5.7)| the source is smooth in e. This is not apparent at a first 
sight because the four last equations have singular sources. But using energy-momentum 
conservation, we see that the projections of the nonphysical divergence (i.e. the conformal 
continuity equations) become : 

(5.9) pj + D'Ji = Q-^op-exp-e^^x^i + ^i^tif + e^oS) 



Ji t + D Sn 



Q^pEi + 2S Ji + Q £ (2S J % - S^) -exJi + efi J { . 



Consequently, inserting |(2.6)| in definition |(2.7)| for t^ uX and using the continuity equa- 
tions, we obtain: 



(5.10) 



(5.11) 



2Qt 



1 



OjO 



D . p - £ ^ p E j + - VSa - -~ 9ij S - e'Q- l E J, + 



-2* Jl 



e 4 ^ 



n- 1 t l s lj --st J ) + -xJ j 



0(ij) 



eD {l J j} + ~g i:j (p tt - Q V^o) + £ 2 pXij + 



s 3 ^ 



(5.12) 
(5.13) 



2Vltij k 



eDyJi] — £ 2 (l~ 1 T l [ i J j ] — e A X[iSj 



1 



D[ip (jj] k + e 2 Vl ( -pg j]k + e 2 (S j]k - -Sg j]k ) ) + 



3 

+e 2 D, ( ,V, 



j]k 



Sgj\k) — & E Q e 2 g k [iJj] — 



e Q, £ g k[i -pg jV + e (Sj]i - -%]«) + e Xk[iJj]- 



Using these expressions it is easy to see that 

~g v n viil = o 
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and 

t n \itj)ni = 0(e 2 ), 
showing that the apparently singular terms are in fact regular. 
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